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ABSTRACT 

Reliable atomic data have been computed for the spectral modeling of the ni- 
trogen K lines, which may lead to useful astrophysical diagnostics. Data sets com- 
prise valence and K-vacancy level energies, wavelengths, Einstein ^-coefficients, 
radiative and Auger widths and K-edge photoionization cross sections. An impor- 
tant issue is the lack of measurements which are usually employed to fine-tune 
calculations so as to attain spectroscopic accuracy. In order to estimate data 
quality, several atomic structure codes are used and extensive comparisons with 
previous theoretical data have been carried out. In the calculation of K pho- 
toabsorption with the Breit-Pauli i?-matrix method, both radiation and Auger 
damping, which cause the smearing of the K edge, are taken into account. This 
work is part of a wider project to compute atomic data in the X-ray regime to 
be included in the database of the popular xstar modeling code. 

Subject headings: atomic processes — atomic data — line formation — X-rays: 
general 
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Introduction 



The improved resolution and sensitivity of current satellite-borne X-ray observatories 
( Chandra and XMM Newton) are allowing the study of previously inaccessible weak spectral 
features of astrophysical interest. In the early stages of these missions, it was realized 
that absorption by near-neutral species was common, and the fact that all charge states 
(except the fully ionized) left identifiable imprints in the X-ray spectrum has proven to be a 
powerful diagnostic. In ner-shell absorption is important in the outflows of Seyfert galaxies 



in terms of both 



( IBehar fc Netzer 



oxygen (jPradhan 



e La ( ISako et al 



2001 



Behar et al. 



200 ll ) and Ka lines of high-Z elements 



2002|), and also of elements in the first row of the periodic table such as 



20001 : 



Behar et al. 



2003 



Garcia et al. 



20051 ) . Furthermore, inner- shell 



absorption of contin uum X-rays from bright galactic sou rces is a useful diagnostic of the 



interstellar medium (lYao et al. 



2009 



Kaastra et al. 



2009|). 



Nitrogen K-shell absorption and emission are detected in X-ray spectra, mostly due to 
the H- and He-like charge states. For i nstance, observations of t he emission lines of N VI 



and N VII in the ejecta of rj Carinae by 



Leutenegger et al. 



(120031 ) have resulted in the lower 



bound N/O > 9 for its nitrogen abundance. This result puts a constraint in the evolution 
of r) Car and is a signature of CNO-cycle processing. The N VI triplet is observed in the 
spectrum of the mag netic B star (3 Cep where it has been used to test magnetically confined 



wind shock models (jFavata et al. 



20091 ). It is found that the plasma is not heated by 



magnetic reconnection and there is no evidence for an optically thick disk at the magnet ic 



equat 



or. Highly ionized emission lines of nitrogen have been observed by 



Miyata et al. 



( 120081 ) on the north-eastern rim of the Cygnus Loop supernova remnant which can be used 



to determine nitrogen abundances, which turn out to be 23% solar. Seyfert 



also can have super-solar N abundances (IBrinkman et al. 



2002 



AravetaL 



galax y outflows 



20071 ). 



Narrow absorption Ka and K/3 lines of N VI have been identified in the warm absorber 
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of the MR 2251- 178 quasar, which po int to a complex velocity field with an outflow of 



ionized material (iRamfrez et al.ll2008l ). K-shell absorption by Li-like N V, which has 
prominent UV lines, or by lower charge states at wavelengths A > 29 A is difficult to 
detect due to the reduced sensitivity of current X-ray instruments towards these longer 
wavelengths. Notable exceptions incl ude : Ka absorption by N V at 29.42 A has been 



reported by ISteenbrugge et al. 



(120051 ) in the outflow of NGC 5548; and N absorption by 



a white dwarf ou tflow has been observed following the outburst of nova V4743 Sagittarii 



( INess et al. 



20031 ). In the latter case, only the H- and He-lines are discussed, but lower 



charge states of N are clearly seen in the spectrum longward of 29 A (see Fig. 3b in 



Ness et al. 



2003L 



The current proliferation of X-ray spectra with high signal-to-noise ratio in astronomical 
archives makes the computation of nitrogen K-shell photoabsorption particularly timely. 



This is a 

data by 



so an additional an d important step in the ong oing effort to compute r eliable atomic 



Garcia et al. 



(120051 ): 



Palmer i et al. 



(12008a 



Witthoeft et al. 


f 


2009 


) for K-line 


deling code (IKallman & Bautista 



200ll ) . Available spectroscopic data for K- vacancy levels of the N isonuclear sequence 
are notably scarce, mainly limited to N V and N VI for which a few levels are listed in 



the NIST d atabase ( 



reported by 



Ralchenko 



Beiersdorfer et al 



et al. 



20081 ) and four measured wavelengths have been 
(119991 ). This shortage of reliable measurements precludes 
empirical corrections to calculated level energies. On the other hand, several relativistic 
methods have been previously used to generate at omic data for n i troge n K-vacancy 



states: the sadd l e-point complex-rotation method (IDavis fc Chung 



Shiu et a 



method (jHsu et al 



2001 



Lin et a 



1991 



2001 



2002 



Yang fc Chung 



rang et al. 



1995 



2005): t 



1989 



Chung 



1990 



Wang fc Gou 



Dirac-Fock ( MCD 



(119841 ): 



Chen! (119861 ): 



approaches such as 



Chen fc Crasemannl (11987 . 



;hose by 



re Ra leigh-Ritz variational 
2006) : and multiconfigurat ion 



H 



1988|); 



at a &: Grant 



Chen et al. 



( 1983 ): 



Hardis et al. 



(jl997l ). Photoabsorption 



- 6- 



cross sections in the near K-edge region of N ions ha ve been obtained by Hartree-Slater 



central-field computations (IReilman &: Manson 



19791 ). where the resonance structure due 



to quasibound states as well as configuration correlations are neglected. The net effect of 
the resonance structure is to fill in the photoionization cross section below the inner-shell 
threshold altering the shape of the K edge. 

In this paper we report on calculations of energy level structure and bound-bound and 
bound-free transition probabilities for the K-shell of nitrogen. The outline of the present 
report is as follows. The numerical methods are briefly described in § 2 while an analysis 
of the results based on comparisons with previous experimental and theoretical values is 
carried out in § 3. The two supplementary electronic tables are explained in § 4 while some 
conclusions are finally discussed in § 5. 



2. Numerical methods 



The numerical approach used here has been fully described in iBautista et al.l (120031 ). 
Level energies, wavelengths, A-co efficients, and radiative and Auger rate s are comp uted 



with the codes autostructure (lEissner et al 



1974J : 



Badnell 



1986 



19971 ) and hfr ((Cowan 



19811 ). For consistency, configuration-interaction (CI) wave functions of the type 



are calculated with the relativistic Breit-Pauli Hamiltonian 



(2) 



where H nr is the usual non-relativistic Hamiltonian. The one-body relativistic operators 

N 



Fib = 22 /"( maSS ) + /n( d ) + /n (SO) 



(3) 



n=l 
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represent the spin-orbit interaction, f n (so), the non- fine-structure mass variation, / n (mass), 
and the one-body Darwin correction, / n (d). The two-body Breit operators are given by 



where the fine-structure terms are g nm (so) (spin-other-orbit and mutual spin-orbit), g nm (ss) 
(spin-spin), and the non-fine-structure counterparts g nm {.css) (spin-spin contact), g nm {d) 
(two-body Darwin), and g nm (oo) (orbit-orbit). It must be pointed out that hfr neglects 
contributions from the two-body term H 2 h of Equation (J4]). 

In HFR, core-relaxation effects (CRE) are always taken into account since each electron 
configuration is represented with its own set of non-orthogonal orbitals optimized by 
minimizing the average configuration energy. In AUTOSTRUCTURE, on the other hand, 
configurations may be represented with either orthogonal or non-orthogonal orbitals which 
then enables estimates of the importance these effects. In the present calculation five 
approximations are considered in order to study the effects of electron correlation, i.e. 
configuration interaction (CI) and CRE, and thus to estimate data accuracy. 

Approximation ASl Atomic data are computed with AUTOSTRUCTURE including CI 
from only the n = 2 complex. CRE are neglected. 

Approximation AS2 Atomic data are computed with AUTOSTRUCTURE including both 
n = 2 CI and CRE. 

Approximation AS3 Atomic data are computed with AUTOSTRUCTURE including CI 
from both the n = 2 and n = 3 complexes. CRE are neglected. 

Approximation HF1 Atomic data are computed with hfr including CI only from the 
n = 2 complex. 
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Approximation HF2 Atomic data are computed with hfr including CI from both the 
n = 2 and n = 3 complexes. 

Photoabsorption cross sections a re obtained wit h the codes Breit-Pauli R -Matrix bprm 



( iBerrington et al. 



1987 



Seaton 



Bar-Shalom et al. 



20011 ) . In bprm, 



19871 ) and hullac 

wave functions for states of an iV-electron target and a colliding electron with total angular 
momentum and parity Jn are expanded in terms of the target eigenf unctions 

Fi(r) 



J IT 



Xi 



J2 c & ■ 



(5) 



The Xi functions are vector coupled products of the target eigenf unctions and the 
angular components of the incident-electron functions; Fi(r) are the radial part of the 
continuum wave functions that describe the motion of the scattered electron; and A is an 
antisymmetrization operator. The functions <&j are bound-type functions of the total system 
co nstructed with target orbit als. The Breit-Pauli re lativistic version has been developed 



by 



Scott & Burke! (119801 ) and 



Scott fe Taylor! (119821 ). but the inclusion of the two-body 



terms (see Equation Sj) is currently in progress, and thus not inclu ded. Auger and radiativ e 



damping are taken into account by means of an optical p otential (IRobicheaux et al. 



Gorczyca fc Badnell 



1996 



Gorczyca fc McLaughlin 



1995 



20001 ) where the resonance energy 



with respect to the threshold acquires an imaginary component. In the present work, 
the iV-electron targets are represented with all the fine structure levels within the n = 2 
complex. It is important to mention that the bprm approach does not allow the inclusion 
of CRE in the photoionization calculations; therefore, both the initial and final states 
correspond to configurations represented with orthogonal orbitals. Thus, the wave functions 
for the target states are those produced with approximation AS1. 



hullac (Hebrew University Lawrence Livermore Atomic Code, 



Bar-Shalom et al. 



20011 ) is a multiconfiguration, relativistic computing pack age that is bas e d on the 



relativistic version of the parametric potential method by 



Klapisch et al. 



(119771 ). and 
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employs a facto rization- interpolation m ethod within the framework of the distorted wave 



approximation (IBar-Shalom et al. 



1988). It includes the Breit interaction for relativistic 



configuration averages and can take into account part of the correlatio n effects by allowing 



differ ent potentials for each group of configurations. Its newest version (IKlapisch fc Busquet 



20091 ). which is used here, incorporates a number of improvements such as important 
corrections to the photoionization subroutines. For the present work we calculated only 
transition energies and direct photoionization cross sections, but hullac can also efficiently 
compute photo-au toionization resonances by means of the isolated resonance approximation 



(jOreg et al. 



199l|), which are subsequently superimposed on the continuum photoionization 
cross section. Moreover, attempts have been made to adapt hullac to calculate the 
quantum interference of resonan ces with the continuum that leads to Fano-type asymmetric 



profiles (IBehar et al. 



2000 . 



20041 ). However, these calculations with hullac are beyond the 



scope of the present paper. 



3. Energy levels 

Energies have been computed for both valence and K-vacancy levels in the five 
approximations delineated in Section [2] and with hullac. A comparison of approximations 
AS1 with AS2 provides an estimate of CRE while those of AS1 with AS3 and HF1 with 
HF2 give measures of out-of-complex CI. Also, level energy accuracy can be bound with a 
comparison of AS2 and HF1, that is, from two physically comparable approximations but 
calculated with two independent numerical codes. 

In Figure dj average energy differences for AS1 vs. AS2 and AS1 vs. AS3 are plotted 
for each ionic species, 3 < N < 7. It may be seen that while core relaxation effects lowers 
the valence-level energies by around 0.5-0.8 eV, it raises by a similar amount those for the 
K-vacancy levels in species with electron number 3 < iV < 5; as a consequence, transition 
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wavelengths for these ions are expected to be shorter due to this effect. Although n = 3 CI 
also causes a lowering (less than 0.7 eV) of the valence level energies (see Figured]), the 
impact on the K- vacancy level energies is more pronounced (as large as 1.5 eV): it mainly 
lowers levels for species with N < 5 and raises those for N > 5. The latter result is mostly 
due to K-vacancy levels from the n = 2 and n = 3 complexes intermixing in the lowly 
ionized members. 

Figure [2] shows the corresponding quantities for HF1 and HF2. As shown in Figure [2l 
the effect of CI on the energies computed with hfr is somewhat different as they are 
decreased (less than 0.5 eV) for both valence and K-vacancy levels, the minimum occurring 
in ions with N = 5 and N = 6. Level energy differences between the AS2 and HF1 data 
sets are within 0.5 eV which is a reliable accuracy ranking of the present level energies. 

Computed level energies are compared with the few spectroscopic measurements 
available (2 < iV < 3) in Table [Q It may be seen that CRE in autostructure 
(approximation AS2) in general reduce differences with experiment. Furthermore, hfr 
and hullac seem to provide better energies than autostructure; differences of hfr 
(hullac) with experiment not being larger than 1.53 eV (1.19 eV). In the He-like system, 
hullac is particularly accurate for triplet states and less so for singlet states. The accurate 
approximation of HF2 is compared with previously computed term energies in Table [2] 
Although we would not quote present term energies with the same number of significant 
figures as previous results, HF2 values are consistently lower. 

Fine-structure level splittings in HF2 can be problematic as shown in Table [3J It 
may be seen that HF2 gives a value in good accord (1%) with beam-foil measurements for 
the Ai?( 5 P 2 , 5 Pi) splitting of the Is2s2p 2 5 P K-vacancy term of N IV but an intolerably 
discrepant one (factor of 2) for AE( 5 P^, 5 P2). On the other hand, both AS2 splittings are in 
reasonable agreement (10%) with experiment and with results obtained with the relativistic 
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Raleigh-Ritz variational method ( 


Yane & Chung 


1995: 


(Hata & Grant 


1983; 


Hardis et al. 


1984 


). This problem 



Hata Sz Grant 



( 119831 ) to be due to the neglect of the Breit interaction which is the case in HF2. 



The most stringent accuracy requirements for the energies come from the capabilities 
of the astronomical instruments which can observe these transitions. Current instruments, 
the Chandra and XMM-Newton gratings, have a resolving power which is nominally e/Ae < 
1000, which imposes a resolution of 0.4 — 0.5 eV in the energy region of the nitrogen K lines 
(400 — 500 eV), essentially the same accuracy we have achieved in the present calculations. 
However, in spectra with good statistical signal-to-noise it is possible to determine line 
centroids to a factor ~3 more accurately than this. Future instruments, principally the 
grating instruments considered for the International X-ray Observatory (IXO), may also 
have resolving power as high as 3000. Knowledge of transition wavelengths and energy 
levels with this precision is needed for truly unambiguous identification of observed features 
and comparison with models. Clearly, precise laboratory measurements are irreplaceable 
requisites in the theoretical fine-tuning of these calculations in order to reduce the current 
uncertainties. 



4. Wavelengths 

The accuracy of computed wavelengths must be determined without a comparison with 
measurements due to the scarcity of the latter. CRE and n = 3 CI in autostructure 
impact wavelengths in an opposite manner to that displayed for the K-vacancy levels in 
Figure [U Specifically, core relaxation on average shortens wavelengths by as much as 
150 mA in ions with 2 < iV < 5 while the effect is less pronounced on the higher members. 
CI increases wavelengths by ~50 mA for N < 4 and decreases them by as much as 150 mA 
for N > 5. CI in hfr in general increases wavelengths with a maximum of 25 mA at iV = 5. 
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Wavelengths computed with AS2 are on average 7 ± 21 rriA greater than HF1. This 
finding can be further appreciated in a comparison with the few available measured 
wavelengths (see Table H|) where computed values are always greater. HF1 differences with 
experiment are the smallest (less than 37 mA) while CRE in autostructure (AS2) also 
reduces discrepancies. 



Wavelengths for K trans i tions i n nitrogen ions have been 



the MCDF method by 



Chenl (|1986h : 



Chen &: Crasemannl (11987 



previo u sly compute d with 



19881 ): 



Chen et al 



(119970 



While reasonable agreement is found with HF2 for species with N = 3 (average difference 
of 13 ± 25 mA) and N = 5 (average difference of —4 ± 33 mA) , puzzling discrepancies are 
fo und for those with N = 4 and N = 6. As shown in Figure [3j the MCDF wavelengths 



of 



Chen fc Crasemannl (119871 ) show the large average difference with respect to HF2 of 



155 ± 45 mA; i.e., on average, they are significantly longer. The inconsistent situation in 
the C-like ion is somewhat different (see Figure H]) where the average difference with HF2 is 
now —45 ± 449 mA showing a very wide and clustered scatter with questionable deviations 
as large as 800 mA. These comparisons lead us to conclude that the wavelengths computed 
with HF2, our best approximation, are accurate to better than 100 mA. 



5. ^-coefficients 

By comparing ^-coefficients computed with approximations AS1 and AS2, CRE 
on the K radiative decay process may be estimated. Discarding transitions subject to 
cancellation which always display large differences, it is found that, for log A > 10, CRE 
generally causes differences not greater than 20%. However, larger discrepancies (54%) 
are found for transitions undergoing multiple electron jumps such as those tabulated in 
Table O ^-coefficients for these peculiar transitions computed with HF1, which should be 
comparable to AS2, are also included in Table [5] finding good agreement (within 15%); 
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in fact, differences between ^-coefficients computed with approximations AS2 and HF1 
are in general within 22%. Furthermore, CI from the n = 3 complex tends to decrease 
^-coefficients with log A > 10, AS3 being on average 16% lower than AS1 and HF2 6% 
lower than HF1. 



MCDF ^-coefficients by 



Chenl (11980 ) 



Chen fc Crasemannl (1987 



1988|); 



ChenetaL 



( 119971 ) agree with HF2 to around 25% except for N = 6 where they are found to be, 



on average, higher by a factor of 4 (see Figure [5]) . 



accuracy of the MCDF ^-coefficients by 



Chen et al- 



i us ou tcome certainly questions the 



( 119971 ) for C-like nitrogen. We find 



that for logA > 10 present ^-coefficients are accurate to within 20% for transitions not 
affected by cancellation. 



6. Radiative widths 

A comparison of AS1 and AS2 radiative widths for the K- vacancy levels indicates 
that CRE are mainly noticeable in the highly ionized species, namely those with N < A, 
where on average the radiative widths are increased by around 10%. On the other hand, 
the inclusion of levels from the n = 3 complex in the CI expansion (AS3) leads to slightly 
reduced radiative widths (less than 7%) with respect to AS1 for the lowly ionized members 
(N > 5). Larger reductions (~20%) are also observed for the lowly ionized species between 
HF2 (which contains CI from the n = 3 complex) and HF1 (which contains CI only from 
the n = 2 complex). 

A remarkable exception is the radiative width of the ls2s 2 2p 3 level in the C-like 
ion. For this level, AS2 gives Aj = ^ Aji = 7.09 x 10 6 s -1 in good agreement with 
HF1 (7.65 x 10 6 s -1 ); however, the AS3 and HF2 radiative widths are respectively 
Aj = 5.31 x 10 8 s -1 and Aj = 3.77 x 10 9 s^ 1 , i.e. around two orders of magnitude larger. 
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The reason for this huge increase when levels from n = 3 complex are included in the CI 
expansion may be appreciated in Table [6J Within the n = 2 complex, the ls2s 2 2p 3 
decays radiatively to the ls 2 2s 2 2p 2 3 Pj ground levels via two spin-forbidden Ka transitions 
which have small ^-coefficients (< 10 7 s _1 ). When the n = 3 complex is taken into account, 
3^2 spin-allowed channels appear which exhibit considerably larger ^-coefficients that 
add up to the quoted enhancement. The observed discrepancy between AS3 and HF2 (a 
factor of 7) are due to severe cancellation in the An ^ transitions. 



Rad iative widths computed for ions with e 



method (IChen 



1986 



Chen Sz Crasemann 



1987 



ectron number 3 < N < 6 with the MCDF 



1988 



Chen et al. 



19971 ) agree with HF2 to 



around 20% except for the C-like species where MCDF is a factor of 4.6 higher. 



7. Auger widths 



Auger widths computed with AS2 and HF1 agree to within 20% but are sensitive to 
both CRE and CI as depicted in Figures [6H3 A comparison of AS1 and AS2 shows that, on 
average, CRE effects increase Auger widths with logA a > 12 linearly as a function of the 
ion electron number N. In autostructure CI also increases Auger widths, particularly 
for highly ionized species (around 28% for ions with 3 < N < 4) while in hfr Auger widths 
are decreased (up to 15% for the lowly ionized members). We believe these contrasting 
outcomes are due to the way orbitals have been optimized in autostructure. 



E xclud in g the C-like ion, M CDF Auger widths by 



(119871 . 



1988 



Chen et al. 



Chenl (1198a ): 



Chen fc Crasemann 



(119971 ) with logA a > 12 in general agree with HF2 to around 20%. 



Larger discrepancies are encountered for a handful of K- vacancy levels listed in Table [7] 
which are mainly caused by level admixture. It may be seen therein that Auger widths 
computed with our different approximations also display a wide scatter thus supporting 
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this diagnostic. For the C-like ion, the MCDF Auger widths are, on average, a factor of 3 
higher than HF2 and thus believed to be of poor quality. 

In Table [HJ Auger widths for the I s2s2p levels of N V co mputed with the Breit-Pauli 



saddle-point complex-rot ation method (IDavis fc Chung 



19891 ) are compared with AS2, HF2, 



and the MCDF values of I Chen! (jl986l ). The agreement with AS2 is within 15% while very 
large differences are found for the HF2 4 P° levels which are most surely due to the neglect 
of the two-body Breit interaction in hfr. The accord with MCDF is around 25% except 
for the Is2s2p 2 P° term where a discrepancy of a factor of 2 is encountered. The latter is 
difficult to explain. 

Auger widths for K- va cancy terms in N IV computed with the Breit-Pauli saddle-point 



complex-rotation method ( 



AS2, HF2, and MCDF flChen fc Crasemann 



Lin et al 



2001 



2002 



Zhang et al. 



20051 ) are compared with 



19871 ) in Table M The level o 



Zhang et al 



agreem ent with 
J2OO5J) for the 



AS2 and HF2 is around 20% except for the values quoted by 
ls2p 3 3 P° and 3 D° terms which are discrepant by about 50%, in contrast with the MCDF 
Auger widths for these two terms which agree with AS2 and HF2 to within 10%. The rest 
of the MCDF Auger widths are in good agreement except for the Is2s2p 2 3 P term which 
has already been singled out in Table [7| as being sensitive to level mixing. 



8. Photoabsorption cross sections 



In Figure [H we show the high-energy photoabsorption cross sections of N I - N V 
computed with the bprm package. Intermediate coupling has been used by implementing 
the AS1 approximation to describe the target wave functions. In order to resolve accurately 
the K-threshold reg;ion, radi a tive a nd Auger damping are taken into account as described 



by 



Gorczyca fc McLaughlin! (120001 ) using the Auger widths calculated with the HF2 
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approximation. For comparison, we have incl uded the photoioniza t ion c ross sections 



Reilman fc Mansonl (119791 ). the latter 



obtained with the hullac code and those by 
calculated in a central-field potential. This comparison shows that the K-threshold energies 
of bprm and hullac are in very good accord (within 1 eV) and the background cross 
sections to within ~ 10%. Note that for the sake of comparison, we have used hullac 
only to compute the direct bo und-free photoioniza t ion cr oss section and not the resonances. 



Background cross sections by 



Reilman &: Mansonl (119791 ) are in excellent agreement with 



bprm for all ions, but K-edge positions and structures are clearly discrepant. The present 
bprm calculations result in smeared K edges due to the dominance of the Auger spectator 
(KLL) channels over the participator ( KLn) channels. This overall behavior is sim ilar to 



that reported in previous calculations (IKallman et al. 



2004 ; 



Witthoeft et al. 



2009 ), and in 



parti cular to that displayed by the corresponding oxygen ions (see Figure 5 in 



Garcia et al. 



20051 ). 



9. Machine-readable tables 

We include two machine readable tables. For nitrogen ions with electron occupancy 
N = 1 — 7 and for both valence and Auger levels, Table [TU] tabulates the spin multiplicity, 
total orbital angular momentum and total angular momentum quantum numbers, 
configuration assignment, energy, and radiative and Auger widths. Fields not applicable 
or not computed are labeled with the identifier "— 9.99E+2". For K transitions, Table fTTl 
tabulates the wavelength, A-coefficient, and gf-value. The printed stub versions list the 
corresponding data only for the Li-like ion (N = 3). The photoionization curves may be 
requested from the corresponding author. 
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10. Discussion and conclusions 

Detailed calculations have been carried out on the atomic properties of K-vacancy 
states in ions of the nitrogen isonuclear sequence. Data sets containing energy levels, 
wavelengths, ^-coefficients and radiative and Auger widths for K-vacancy levels have 
been computed with the atomic structure codes hfr and autostructure. High-energy 
photoionization and photoabsorption cross sections for members with electron occupancies 
N > 3 have been calculated with the bprm and hullac codes. 

Our best approximation (HF2) takes into account both core-relaxation effects and 
configuration interaction from the n = 3 complex. By comparing results from different 
approximations with previous theoretical work and the few spectroscopic measurements 
available, we conclude that level energies and wavelengths for all the nitrogen ions 
considered in the present calculations can be quoted to be accurate to within 0.5 eV and 
100 niA, respectively. The accuracy of ^-coefficients and radiative and Auger widths is 
estimated at approximately 20% for transitions neither affected by cancelation nor strong 
level admixture. Outstanding discrepancies are found with some MCDF data, in particular 
wavelen gths and ^-coeffici ents for the C-like ion which are believed to be due to numerical 



error by 



Chen et al. 



(119970 . 



We have also presented detailed photoabsorption cross sections of nitrogen ions in 
the near K-threshold region. Due to the lack of previous experimental and theoretical 
results, we have also performed simpler calculations using the hullac code in order to 
check for consistency and to estimate accuracy. Comparison of bprm and hullac indicates 
that present K-threshold energies are accurate to within 1 eV. However, background cross 



Reilman fc Mansonl (119791 ) in a 



sections are in better agreement with those computed by 
central-field potential for all ions except the Li-like system. These photoabsorption cross 
sections and their structures are similar to those displayed by ions in the same isoelectronic 
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2005|). 



seque nce (jWitthoeft et al.ll2009l ). especially to the corresponding oxygen ions (IGarcia et al. 



The present atomic data sets are available on request and will be incorporated in the 
xstar modeling code in order to generate improved opacities in the nitrogen K-edge region, 
which will lead to useful astrophysical diagnostics such as those mentioned in § 1. 
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3 4 5 6 

Electron Number 



Fig. 1. — Average level energy differences with respect to approximation AS1 for ions of 
the nitrogen isonuclear sequence with electron number 3 < N < 7. Filled circles: valence 
level energies computed with AS2. Open circles: valence level energies computed with AS3. 
Filled squares: K- vacancy level energies computed with AS2. Open squares: K- vacancy level 
energies computed with AS3. 
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4 5 6 

Electron Number 



7 



Fig. 2. — Average level energy differences with respect to approximation HF1 for ions of the 
nitrogen isonuclear sequence with electron number 3 < iV < 7. Open circles: valence level 
energies computed with HF2. Open squares: K-vacancy level energies computed with HF2. 
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Fig. 3. — Wavelength differences between MCDF (I Chen fc Crasemann 
N IV. An average difference of 155 ± 45 mA is observed. 



19871 ) and HF2 for 
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Fig. 4. — Wavelength differences between MCDF (jChen et al. 
Differences as large as 800 mA are observed. 
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1222) and HF2 for N II. 
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Fig. 5. — A-coefficient differences (s 1 ) between MCDF ( jChen et al 
It is found that MCDF is on average higher by a factor of 4. 



19971 ) and HF2 for N II. 
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Fig. 6. — Average percentage difference between Auger widths (log A a > 12) computed with 
the AS2 and AS1 approximations (filled squares) and with AS3 and AS1 (open squares). 
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Fig. 7. — Average percentage difference between Auger widths (\ogA a > 12) computed with 
the HF2 and HF1. 
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Fig. 8. — High-energy photoabsorption cross sections for ni trogen ions in the K-edge r egion. 
Solid curve: bprm. Dotted curve: hullac. Dashed curve: iReilman fc Mansonl (I1979I ). 
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Table 1. Measured and computed K- vacancy level energies (eV) 



N 


Level 


E 






E(Th) - 


£7 (Expt) 




Expt 


AS1 


AS2 


AS3 


HFl HF2 


HULLAC 


2 


ls2s 3 Si 


419.80 


-1.76 


-1.34 


-2.35 


-0.59 


0.01 


2 


ls2p 3 P° 


426.30 


-1.85 


-0.93 


-2.36 


-0.56 


0.09 


2 


ls2p 3 P° 


426.30 


-1.85 


-0.93 


-2.36 


-0.54 


0.11 


2 


ls2p 3 P° 


426.33 


-1.86 


-0.94 


-2.37 


-0.53 


0.13 


2 


ls2s % 


426.42 


-1.61 


-0.21 


-2.15 


0.03 


0.69 


2 


ls2p X P° 


430.70 


-1.64 


-0.85 


-2.24 


-0.24 


0.67 


3 


ls( 2 S)2s2p( 3 P°) 4 P° /2 


414.61 


-2.49 


-1.97 


-3.64 


-1.48 -1.53 


0.42 


3 


ls( 2 S)2s2p( 3 P°) 2 P° /2 


421. 52 a 


-1.72 


-0.80 


-2.18 


-0.51 -0.62 


0.83 


3 


ls( 2 S)2p 2 ( 3 P) 4 P 5/2 


425.70 


-2.54 


-1.64 


-3.37 


-1.31 -1.33 


1.19 



Note. — Experimental level e nergies (relative to the ion ground state) from the NIST 
database (IRalchenko et al.ll2008l ). 

a Derived from the wavelength measurement of iBeiersdorfer et all (119991 ). 
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Table 2. Computed total energies (au) for .fT-vacancy terms 



N Term HF2 Other theory 



3 


Is2s2p 4 P° 


-33.240 


-33.192008 a , 


-33.192204 b 


3 


ls( 2 S)2s2p( 3 P°) 2 P° 


-32.952 


-32.919222 a 




3 


ls( 2 S)2s2p( 1 P°) 2 P° 


—32.777 


—32 768550* 




4 


ls2s 2 2p 3 P° 


-36.219 


-36.171232 c , 


-36.173064 d 


4 


Is2s2p 2 3 S 


-35.646 


-35.615357 d 




1 


Is2s2p 2 3 P 


-35.842 


-35.788866 d 




4 


Is2s2p 2 3 P 


-35.542 


-35.536868 d 




4 


Is2s2p 2 3 D 


-35.815 


-35.785042 d 




4 


ls2s 2 2p ip 


-36.081 


-36.036967° 




4 


ls2p 3 1 P° 


-35.086 


-35.082958° 




4 


Is2s2p 2 1 D 


-35.594 


-35.583448 f 




4 


Is2s2p 2 1 P 


-35.445 


-35.435017 f 




4 


Is2s2p 2 X S 


-35.425 


-35.415017 f 




4 


Is2s2p 2 5 P 


-36.160 


-36.0934586 s 


, -36.0934407 


4 


ls2p 3 5 S° 


-35.596 


-35.5414665 s 


, -35.5413248 


4 


ls2p 3 3 P° 


-35.214 


-35.204065 1 




1 


ls2p 3 3 D° 


-35.387 


-35.36660P 





Breit— Pauli saddle-point complex-rotation method ( Davis fc Chung 



1989) 



b Relativistic Raleigh-Ritz variational method I Hsu et al 



1991) 



c Breit-Pauli saddle-point complex-rotation method (|Chune 



1990) 



d Breit-Pauli saddle-point complex-rotation method ( Lin et al 



2001) 



e Breit-Pauli saddle-point complex-rotation method (ILin et al 



2002) 



f Breit-Pauli saddle-point complex-rotation method ( Shiu et al 



2001) 



BRelativistic Raleigh-Ritz variational method ( Wang fc Gou 



2006) 



h Relativistic Raleigh-Ritz variational method I Yang &: Chung 



1995) 



'Brcit-Pauli saddle-point complex-rotation method rtZhang et al 



2005) 
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Table 3. Level splittings (cm 1 ) for the Is2s2p 2 5 P K-vacancy term of N IV 



Level splitting Expt a HF2 AS2 Other theory 



A£( 5 P 2 , 5 Pi) 127± 1 126 119 127.07\ 126.9 C , 129.19 d , 128.94 s 
A£( 5 P 3 , 5 P 2 ) 79.5± 0.8 188 70 78.49 b , 78.45 c , 86.72 d , 86.58 e 



'Beam-foil measurements by iBerry et al.l (119821 ) 



b Relativistic Raleigh-Ritz variational method (IWang &: Gou 



20061 ) 



BRelativistic Raleigh-Ritz variational method (jYang & Chung! Il995l ) 



d MCDF-EAL calculation of 



Hata fc Grant 



(1983) 



°MCDF calculation of 



Hardis et al. 



(11984h 
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Table 4. Experimental and theoretical wavelengths (A) for nitrogen ions 



N 


Lower level 


Upper level 


A(Expt a ) 




A(Th) - A(Expt) 












HF1 


AS1 


AS2 


AS3 


HULLAC 


2 


Is 2 % 


ls2s 3 Si 


29.5321(26) 




0.1266 


0.0971 


0.1683 


0.0012 


2 


Is 2 % 


ls2p 3 P° 


29.0835(26) 


0.0374 


0.1274 


0.0642 


0.1627 


-0.0074 


2 


Is 2 % 


ls2p X P° 


28.7861(22) 


0.0167 


0.1108 


0.0576 


0.1515 


-0.0444 


3 


ls 2 2s 2 S 1/2 


Is2s2p 2 P° /2 


29.4135(37) 


0.0353 


0.1208 


0.0560 


0.1550 


-0.0573 



a Wavelength measurements by 



Beiersdorfer et al 



(1999) 



Table 5. Core relaxation effects on K radiative decay 



N 


j 


i 






A 3i (s- 1 ) 




AS1 


AS2 


HF1 


3 


ls2s 2 2 S 1/2 


ls 2 2p 2 


po 
r l/2 


3.85E+10 


5.36E+10 


4.63E+10 


3 


ls2s 2 2 S 1/2 


ls 2 2p 2 


po 

r 3/2 


7.62E+10 


1.06E+11 


9.18E+10 


4 


ls2s 2 2p X P° 


ls 2 2p 2 


X D 2 


9.72E+10 


1.36E+11 


1.24E+11 


5 


ls2s 2 2p 2 2 P 1/2 


ls 2 2p 3 


2po 
r l/2 


1.17E+10 


1.80E+10 


1.78E+10 


5 


ls2s 2 2p 2 2 P 3/2 


ls 2 2p 3 


2po 
r 3/2 


1.54E+10 


2.36E+10 


2.32E+10 
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Table 6. Radiative decay routes of ls2s 2 2p 3 5 S£ in N II 



3 i 


Aji 


(s- 1 ) 




AS3 


HF2 


ls2s 2 2p 3 5 S° ls 2 2s 2 2p 2 3 P X 


2.04E+6 


8.57E+6 


ls 2 2s 2 2p 2 3 P 2 


2.97E+6 


2.53E+7 


ls 2 2s2p 2 3d 5 Pi 


5.07E+7 


3.43E+8 


ls 2 2s2p 2 3d 5 P 2 


8.36E+7 


5.70E+8 


ls 2 2s2p 2 3d 5 P 3 


1.17E+8 


7.97E+8 


ls 2 2p 3 3p 5 Pi 


7.85E+7 


6.99E+8 


ls 2 2p 3 3p 5 P 2 


1.31E+7 


1.16E+9 


ls 2 2p 3 3p 5 P 3 


1.83E+8 


1.63E+8 



J2i A ji 



5.31E+8 3.77E+9 
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Table 7. Discrepant Auger rates (s x ) 



N 


Level 


ASl 


AS2 




AS3 


HFl 


HF2 




MCDF a 


3 


ls( 2 S)2s2p( 3 P°) 2 P° /2 


1.51E+13 


5.81E4 


-12 


2.15E4-13 


8.61E4-12 


7.63E4 


12 


1.41E4-13 


3 


ls( 2 S)2s2p( 3 P°) 2 P° /2 


1.43E+13 


5.27E4 


-12 


2.02E4-13 


8.29E4-12 


7.30E4 


12 


1.35E4-13 


4 


ls( 2 S)2s2p 2 ( 4 P) 3 P 


3.09E+13 


1.99E4 


-13 


4.41E4-13 


1.86E4-13 


1.77E4 


13 


3.73E4-13 


4 


ls( 2 S)2s2p 2 ( 4 P) 3 Pi 


3.08E+13 


1.98E4 


-13 


4.40E4-13 


1.86E4-13 


1.76E4 


13 


3.68E4-13 


1 


ls( 2 S)2s2p 2 ( 4 P) 3 P 2 


3.04E+13 


1.95E4 


-13 


4.36E4-13 


1.85E4-13 


1.75E4 


13 


3.55E4-13 


4 


ls( 2 S)2s2p 2 ( 2 P) 1 P 1 


2.26E+13 


2.88E4 


-13 


2.62E4-13 


2.46E4-13 


2.34E4 


13 


1.43E4-14 


1 


ls( 2 S)2s2p 2 ( 2 S) 1 Sq 


1.34E+14 


1.33E4 


-14 


1.77E4-14 


1.28E4-14 


1.19E4 


14 


1.77E4-13 


5 


ls( 2 S)2s2p 3 ( 5 S°) 4 S° /2 


3.90E+13 


2.98E4 


-13 


4.10E4-13 


2.54E4-13 


2.48E4 


13 


3.89E4-13 


5 


ls( 2 S)2s2p 3 ( 3 S°) 2 S° /2 


2.34E+13 


3.09E4-13 


2.88E4-13 


5.30E4-13 


1.31E4-14 


3.49E4-13 


5 


ls( 2 S)2s2p 3 ( 1 P°) 2 P° /2 


1.37E+14 


1.67E4 


-14 


1.56E4-14 


1.27E4-14 


3.48E4 


13 


1.52E4-14 



MCDF computations bv lChenl \l98Q ): IChen fc Crasemannl p987l . Il988h 



Table 8. Auger energy widths (au) for Is2s2p levels in N V 



Level 




AS2 




HF2 




Other theory 




ls( 2 S)2s2p( 3 P°) 


4po 
r l/2 


1.46E- 


-08 


3.48E- 


-10 


1.50E-08 a , 1.532E 


-08 b 


ls( 2 S)2s2p( 3 P°) 


4po 
r 3/2 


4.53E- 


-09 


8.83E- 


-10 


4.98E-09 a , 3.952E 


-09 b 


ls( 2 S)2s2p( 3 P°) 


4po 
r 5/2 


4.26E- 


-10 






4.587E-10 b 




ls( 2 S)2s2p( 3 P°) 


2po 


1.32E- 


-04 


1.79E- 


-04 


3.31E-04 a , 1.54E- 


04 b 


ls( 2 S)2s2p( 1 P°) 


2po 


1.68E- 


-03 


1.47E- 


-03 


1.14E-03 a , 1.53E- 


03 b 



"MCDF calculations by (IChen 



1986|) 



b Breit-Pauli saddle-point complex-rotation method (IDavis &: Chung 



1989) 
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Table 9. Auger widths (meV) for K- vacancy terms in N IV 



Term AS2 HF2 Other theory 



ls2s 2 2p 3 P° 


84.1 


71.8 


77.4 a . 


79.0 b 


ls( 2 S)2s2p 2 ( 4 P) 3 P 


12.9 


11.6 


23.8 a . 


10.8 b 


ls( 2 S)2s2p 2 ( 2 D) 3 D 


66.9 


62.6 


53.4 a 


57. 7 b 


ls( 2 S)2s2p 2 ( 2 S) 3 S 


31.2 


29.7 


26. 7 a 


28.6 b 


ls( 2 S)2s2p 2 ( 2 P) 3 P 


66.4 


55.5 


58.5 a 


55.0 b 


ls2s 2 2p 1 P° 


55.2 


48.3 


53.8 a 


58 c 


ls2p 3 1 P° 


44.9 


47.0 


43.8 a 


43 c 


ls2p 3 3 D° 


80.7 


80.3 


75.0 a 


53.87 d 


ls2p 3 3 P° 


47.4 


47.5 


45.5 a 


34.16 d 



MCDF method flChen k Crasemannlll987h 



5 Breit-Pauli sadd 



method (ILin et al 



2001) 



e-point complex-rotation 



'Breit-Pauli sadd 



method (ILin et al. 



20021 ) 



e-point complex-rotation 



d Breit -Pauli saddle-p oint complex-rotation 



method (jZhang et al. 



2005|) 
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Table 10. Valence and Auger levels for nitrogen ions 



N 


i 


25 + 1 


L 


2J 


Configuration 


Energy 
eV 


At 
s- 1 


Aa. 
s- 1 


Q 
O 


i 
i 


O 

z 


n 
U 


i 

1 


1 nOOn TCI /O 

iszzs Zoi/z 


n nnnn 
U.UUUU 


Ci OOT7 i no 


Cl OOT7 i no 

— y.yyuj-l-uz 


3 


2 


2 


1 


1 


lsZzp z.roi/2 


n n a £n 
9. 9459 


n qott i no 


— 9.99rij+02 


3 


3 


2 


1 


3 


ISZZp Z.rOo/2 


n n r 7 r 7/i 
9.9/ /4 


q q XT i no 


n nnT7 i no 

— y.yyUj+uz 


3 


A 

4 


2 





1 


ISZSz Zol/Z 


A 1 f\ 1 KiiO 


1 OOTT 1 1 1 

l.ZZrL+11 


O 7(117 1 1 Q 

9. (9-k+lo 


3 


5 


4 


1 


1 


ls(2S)2s2p(3Po) 4Pol/2 


413.0288 


5.05E+06 


1.44E+07 


3 


6 


4 


1 


3 


ls(2S)2s2p(3Po) 4Po3/2 


413.0471 


1.27E+07 


3.65E+07 


3 


7 


4 


1 


5 


ls(2S)2s2p(3Po) 4Po5/2 


413.0778 


0.00E+00 


1.76E+07 


3 


8 


2 


1 


1 


ls(2S)2s2p(3Po) 2Pol/2 


420.8755 


1.65E+12 


7.63E+12 


3 


9 


2 


1 


3 


ls(2S)2s2p(3Po) 2Po3/2 


420.8968 


1.65E+12 


7.30E+12 


3 


10 


4 


1 


1 


ls(2S)2p2(3P) 4P1/2 


424.3267 


7.99E+08 


5.54E+06 


3 


11 


4 


1 


3 


ls(2S)2p2(3P) 4P3/2 


424.3446 


8.06E+08 


4.65E+08 


3 


12 


4 


1 


5 


ls(2S)2p2(3P) 4P5/2 


424.3741 


8.15E+08 


2.85E+09 


3 


13 


2 


1 


1 


ls(2S)2s2p(lPo) 2Pol/2 


425.6391 


1.63E+11 


6.06E+13 


3 


14 


2 


1 


3 


ls(2S)2s2p(lPo) 2Po3/2 


425.6545 


1.55E+11 


6.10E+13 


3 


15 


2 


2 


3 


ls(2S)2p2(lD) 2D3/2 


429.2436 


8.33E+11 


1.07E+14 


3 


16 


2 


2 


5 


ls(2S)2p2(lD) 2D5/2 


429.2442 


8.32E+11 


1.07E+14 


3 


17 


2 


1 


1 


ls(2S)2p2(3P) 2P1/2 


430.3236 


2.68E+12 


1.20E+08 


3 


18 


2 


1 


3 


ls(2S)2p2(3P) 2P3/2 


430.3597 


2.68E+12 


4.50E+10 


3 


19 


2 





1 


ls(2S)2p2(lS) 2S1/2 


437.1874 


7.86E+11 


1.60E+13 



Note. - 



A complete machine-readable version of this table may be accessed on-line. 
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Table 11. Radiative K-transition data nitrogen ions 



N 


3 


i 


Wavelength 


A-coefficient 


g/-value 








0.1 nm 


s- 1 




3 


4 


2 


on d'vdo 

oU.979o 


4.Uo-h+lU 


1 1 r ~rT?> no 


3 


4 


3 


on nooo 


no tt 1 1 in 


dot? no 


3 


5 


1 


on moo 

30.0183 


5.05-b+Ob 


l.dbhj— 0b 


3 


6 


1 


on n i 7n 

30.0170 


1 0717 1 nT 

1.27t+U7 


b.87hj— 0b 


3 


8 


1 


on a et of 

29.458b 


l.D0il(+12 


a onTT 1 ni 


3 


9 


1 


29.4571 


1 r TTi 1 1 

l.DoliiH-lz 


c 1 TT 1 ni 


3 


10 


2 


on <\ n 

29.9203 


9.70-h+Ob 


Z. Dili/ — (Jb 


3 


10 


3 


on < \ c 

29.922b 


2.3b-b+U5 


c /i T? no 


3 


11 


2 


29.9191 


t ni TT' \ r\ a 

7.01t+04 


t'ttt 1 no 

3.77tj— 08 


3 


11 


3 


29.9213 


1.58E+07 


8.49E-0b 


3 


12 


3 


29.9192 


2.14E+07 


1.72E-05 


3 


13 


1 


29.1289 


1.62E+11 


4.12E-02 


3 


14 


1 


29.1279 


1.54E+11 


7.83E-02 


3 


15 


2 


29.5695 


7.16E+11 


3.76E-01 


3 


15 


3 


29.5717 


1.17E+11 


6.12E-02 


3 


16 


3 


29.5717 


8.32E+11 


6.55E-01 


3 


17 


2 


29.4935 


1.79E+12 


4.68E-01 


3 


17 


3 


29.4957 


8.89E+11 


2.32E-01 


3 


18 


2 


29.4910 


4.24E+H 


2.21E-01 


3 


18 


3 


29.4932 


2.26E+12 


1.18E+00 


3 


19 


2 


29.0197 


2.56E+11 


6.46E-02 


3 


19 


3 


29.0218 


5.28E+11 


1.33E-01 



Note. — A complete machine-readable version of this tabic 
may be accessed on-line. 



